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Rhythm Estimator (Rytmeestimator) . ^^•^ 



This invention relates to an apparatus for proc^sing 
auscultation signals, comprising a bias processor for re- 
ceiving an auscultation signal and providing a biased 
auscultation signal; said bias processor comprising an 
envelope detector; an estimator for calculating a signal 
representative of an estimated rhythm of the auscultation 
signal . 



In recent years physicians have had an impressive arsenal 
of instrumentation at their disposal for the diagnosis of 
cardiovascular diseases. One such instrument is the well 
known stethoscope used to detect sounds originating from 

15 the heart and adjacent large vessels. With the introduc- 
tion of the stethoscope, auscultation grew into an art 
providing the clinician with valuable information of the 
functional integrity of the heart. Nowadays it is possi- 
ble to process the information residing in the ausculta- 

20 tion signal electronically by use of embedded knowledge 
obtained by clinical research. Therefore it is possible 
to relieve the clinician from trivial tasks. One such 
task may be to estimate the heart rate of a beating heart 
based on the auscultation signal. 

25 

The auscultation signal comprises characteristic sounds 
that originate from the contraction and expansion of the 
heart in . the process of pumping blood, i.e. sounds cause*d 
by vibrations from acceleration and deceleration of 

30 blood. When a beating heart is auscultated it is not the 
beat frequency that is heard, but rather the characteris- 
tic sounds- Further, for mid-aged healthy persons in nor- 
mal conditions, the characteristic sounds denoted SI and 
S2 related to a contraction and expansion of the heart, 

35 respectively, are separated in time by a sufficiently 
large interval that allows for perception of the charac- 
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teristic sounds as individual sounds, i.e. the character- 
istic sounds, are typically separated by more than 150- 
200 milliseconds for adults and about 50 milliseconds for 
children. Typically, the characteristic sounds have a 
5 continuous frequency spectrum distributed between ap- 
proximately 15 Hz and 100 Hz and peaking at 30-40Hz. 

However, the characteristic sounds may be superimposed 
with sounds originating from different types of vibra- 

10 tions in the heart. These sounds are the so-called mur- 
murs, that may be caused by mechanisms such as blood tur- 
bulence or flutter of structures. Murmurs may be divided 
into two primary groups: a first group comprise random 
^noisy' sounds and a second group comprise ^musical' 

15 sounds with a continuous narrow-band frequency spectrum. 
Experience has shown that in a human being the ^musical' 
sounds are typically related to the heart beat frequency. 
The frequency spectrum of murmurs may be distributed up 
to 500-2000HZ. 

20 

The characteristic sounds may also be disturbed by spuri- 
ous contractions of the heart appearing at certain or 
random points in time. Further, the characteristic sounds 
may also appear at random points in time. This phenomenon 
25 is known as arrhythmia. 

When the sound of a beating heart including murmurs, if 
any, is auscultated by means of a transducer placed on 
e.g. the chest of a patient, then low-frequent rumble 
30 noise will be picked up by the transducer. This rumble 
noise may have a frequency content primarily distributed 
from just below the frequency of the heart rate towards 
DC, 

35 In severe cases these disturbances are able to destroy 
the characteristic sounds in the auscultation signal com- 
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pletely or at least disturb the auscultation signal to 
such a degree that it may be impossible to estimate even 
an approximated heart rate by use of the prior art tech- 
nology, 

5 

However, the human perception of the auscultation signal 
is very good at diminishing the disturbances or artefacts 
and pay attention to the features in the auscultation 
signal that provides information for estimating the heart 
10 rate. Thereby it is superior to the prior art technology 
in the above mentioned severe cases. 

This is because the human hearing system enhances high 
frequencies relative to low frequencies in the frequency 

15 range of interest i.e, below approximately 2000Hz. The 
human perception of sound amplitudes may be modelled by 
means of A-weighing. A-weighing may be approximated as a 
second order band-pass filter with a centre frequency of 
2000Hz*. Thereby it reduces the low frequency rumble noise 

20 and enhances the high frequency components providing well 
defined information of the relative location of the char- 
acteristic sounds with respect to time. This well defined 
information in the form of characteristic heart sounds is 
easily recognizable from repetition to repetition and 

25 therefore easy to count. However, the value of the human 
hearing system has not been taken into account in the 
signal processing used in systems for estimating the 
heart rate automatically. This may be due to the preju- 
dice that enhancement of the relative high frequencies in 

30 the human hearing system also enhances noise and thereby 
decreases the signal-to-noise ratio. 

US patent 4,972,841 discloses an electronic stethoscope 
with a pulse rate display for displaying an estimated 
35 pulse rate. In the stethoscope a transducer converts so- 
called Korotkoff sounds into electrical signals. 
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Korotkoff sounds are sounds created* by vibrating walls of 
collapsed vessels as blood is just allowed to flow 
through the vessel. The electrical signals are amplified 
and fed to a counter in which the detected pulse rate per 
5 unit time is calculated and the result is displayed as a 
digital pulse rate. 

However, this patent is not concerned with estimation of 
heart rates originating from humans where the measured 
10 sound may be very different from expected normal heart 
sounds. The patent does not disclose any means adapted to 
suppress artefacts that may be confused with or mixed up 
with the information in the auscultation signal relating 
to the heart rate. 

15 

US patent 4,436,096 discloses a digital heart rate 
stethoscope. A transducer produces electrical signals 
representing the heart sound of a human patient and a 
filter is connected to receive the electrical signals. 

20 The filter is a narrow-band band-pass filter whose pass- 
band is centered on a characteristic heart sound fre- 
quency e.g. 33 Hz. The filter is thereby adapted to im- 
prove the signal-to-noise ratio and enables the trans- 
ducer to be used over a patient's clothing. The filtered 

25 signal is converted into pulses which are counted in a 
counter and displayed as a digital heart rate indication. 

Again, this patent is not concerned with estimation of 
heart rates originating from humans where the measured 
30 sound may be very different from the expected normal 
heart .sounds . Neither does this patent disclose any means 
adapted to suppress artefacts that may be confused with 
or mixed up with the information in the auscultation sig- 
nal relating to the heart rate. 



35 



5 



International Application WO 97/00045 discloses a stetho- 
scope for examining the characteristic of sound waves 
made by a beating heart. The stethoscope comprises means 
for establishing a threshold wave form beyond which cy- 
5 clic reference complexes appear and are identified, and a 
time measurement device for measuring the time elapsed 
between similar moments on adjacent reference complexes 
as an indicator of rate of heart beat, 

10 In this stethoscope the heart rate is estimated based on 
the assumption that there is substantially no noise be- 
tween the heart sounds originating from contraction and 
expansion of the heart. The stethoscope is therefore not 
suited for use of estimating heart rates of persons suf- 

15 fering from heart murmur or arrhythmia. 

Other apparatuses for estimating the heart rate use means 
for determining the envelope of the auscultation signal. 
The envelope is determined by taking the absolute value 
20 of the auscultation signal and band-pass filtering the 
absolute value signal . 

From the prior art it is clear that there is not dis- 
closed any method or apparatus that uses the knowledge of 

25 characteristic features for typical auscultation signals 
where murmurs and arrhythmia may be present', in order to 
obtain a more precise and robust heart rate estimation. 
Further, the strength of the bias provided by the human 
hearing system in the perception of auscultation signals 

30 is not used in the prior art in an apparatus for reliev- 
ing the clinician of the trivial task of calculating the 
heart rate. 

The primary object of the invention is to provide an ap- 
35 paratus for processing auscultation signals in order to 
extract information representative of the rhythm of a 
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beating heart such that random noise and spontaneous ar- 
tefacts are suppressed. 

This is achieved, when the apparatus is characterised in 
5 that the estimator is adapted for selecting a part of the 
biased auscultation signal and measuring the conformity 
between the part of the biased auscultation signal and 
the entire auscultation signal. 

10 Consequently, it is possible to suppress artefacts in the 
auscultation signal that may be confused or mixed up with 
the information relating to the heart rate due to an 
overlap of the frequency spectrum of the artefacts and 
the frequency spectrum of the heart rate information. 

15 There is thereby obtained a signal in which the heart 
rate is very well defined. 

In an expedient apparatus for processing auscultation 
signals, the estimator calculates a cross-correlation 
20 • function. This cross-correlation function supplies a 
well-defined signal having large signal values at good 
conformity and small signal values at poor conformity. 

When the part of the biased auscultation signal repre- 
25 sents one of a succession of cycles of the biased auscul- 
tation signal, the part may represent a reference cycle 
that is found ideal in some aspects. 

The apparatus may further comprise a post processor 
30 adapted for locating extreme values in the signal repre- 
senting the conformity of the biased auscultation signal. 
Hereby, the points in time in which the extreme values 
are located may be a direct measure of the rhythm. 

35 When the post processor is further adapted for validating 
the quality of the signal representing the conformity of 
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the biased auscultation signal, then the quality of the 
estimated rhythm may be taken into consideration when the 
estimated rhythm rate is interpreted. 

5 The quality may be validated by verifying at least one of 
the following three items: 

a) time differences between located extreme values must 
be within predetermined limits; 

b) minimum and maximum time differences in proportion to 
10 the mean of the time differences must be within prede- 
termined limits; 

c) the magnitude of the result of the correlation at the 
extreme values location must be within predetermined 
limits . 

15 Then a preferred quality level may be defined by means of 
the preferred limits that are related directly to the 
statistics of rhythm estimation. 

The bias processor may comprise a filter for calculating 
20 an A-weighted version of the auscultation signal or "an 
approximated A-weighted version of the auscultation sig- 
nal. Simulations have shown that an A-weighted version of 
the auscultation signal comprises sufficient information 
to estimate the heart rate in a consistent way in a vari- 
25 .ety of situations. 

When the A-weighted version of the auscultation signal is 
calculated by means of an approximation corresponding to 
a double differentiation of the auscultation signal in a 
30 predetermined frequency range, the filtering process re- 
quires a minimum of computational effort. 

The bias processor may comprise an envelope detector for 
calculating a first envelope signal representing the en- 
35 velope of the A-weighted version of the auscultation sig- 
nal. This is expedient because the envelope signal of the 



A-weighted version of the auscultation signal comprises 
information that is very feasible for estimating the 
heart rate even when murmurs are present in the ausculta- 
tion signal. 

5 

When the bias processor comprises an envelope detector 
for calculating a second envelope signal representing the 
envelope of the auscultation signal, the envelope signal 
of the auscultation signal is very feasible for estimat- 
10 ing the heart rate when murmurs are not present in the 
auscultation signal. 

The bias processor may comprise an adaptive band-pass 
filter for filtering the envelope signals; said band-pass 

15 filter at least having an upper and a lower pass-band, 
respectively, to select from; said adaptive band-pass 
filter comprising a controller selecting the lower pass- 
band when a relatively large fraction of a signal input 
to the band-pass filter is low-frequent, and selecting 

20 the upper pass-band when a relatively low fraction of a 
signal input to the band-pass filter is low-frequent* 
Thus, it is possible to estimate the rhythm- in the aus- 
cultation signal over a relatively wide range of frequen- 
cies while it is avoided that the signal-to-noise ratio 

25 is ruined, i.e. it is possible preserve a relatively nar- 
row pass-band of the band-pass filter. 

The bias processor may comprise means for combining the 
envelope signals or the band-pass filtered envelope sig- 
30 nals; said means comprising auto-scaling means and summa- 
tion means providing the biased auscultation signal. It 
is thereby possible to combine differently biased signals 
in order to estimate the heart rate of a broader spectrum 
of types of auscultation signals. 
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The auscultation signal may comprise samples that arrive 
at a sample rate, and the apparatus may comprise a syn- 
chronous processor operating at a rate corresponding to 
the sample rate, and further comprising an asynchronous 
5 processor operating at time intervals that are initiated 
by a request. Hereby, it is possible to acquire new sam- 
ples of the auscultation signal in real time while a cor- 
relation function operating on prior sample values of the 
auscultation signal is calculated. 

10 

A second object of the invention is to provide an appara- 
tus for processing auscultation signals in order to ex- 
tract information relating to the repeated interdependent 
contraction and expansion of a beating heart, such that 
15 the information in the characteristic sounds and in mur- 
murs related to the heart rate are enhanced. 

This is achieved, when the apparatus is characterised in 
that the bias processor comprises a filter having a fre- 
20 quency response corresponding to an A-weighing or an ap- 
proximated A-weighing, at least for a frequency range of 
interest. 

Consequently, it is possible to estimate the heart rate 
25 of a beating heart, based on the characteristic sounds of 
the beating heart and murmurs related to the heart rate. 
Thus, the heart rate may be estimated based on murmurs 
relating to the heart rate if the characteristic sounds 
do not provide sufficient information for estimating the 
30 heart rate. 

When the auscultation signal is filtered with a filter 
having a frequency response corresponding to a double 
differentiation, the computational effort required to 
35 calculate the A-weighted version of the auscultation is 
reduced . 
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The frequency range of interest may span the frequency 
range below 2000Hz. The characteristic heart sounds and 
murmurs related to the heart rate comprise frequency com- 
5 ponents located in this frequency range and provide suf- 
ficient information in this range for estimating the 
rhythm of a beating heart. 

In an expedient apparatus for processing auscultation 
10 signals, the estimator is adapted for selecting a part of 
the biased auscultation signal and measuring the confor- 
mity between the part of the biased auscultation signal 
and the entire auscultation signal. 

15 In an expedient apparatus for processing auscultation 
signals, the estimator calculates a cross-correlation 
function. . This cross-correlation function supplies a 
well-defined signal having large signal values at good 
conformity and small signal values at poor conformity, 

20 

When the part of the biased auscultation signal repre- 
sents one of a succession of cycles of the biased auscul- 
tation signal, the part may represent a reference cycle 
that is found ideal in some aspects. 

25 

The apparatus may further comprise a post processor 
adapted for locating extreme values in the signal repre- 
senting the conformity of the biased auscultation signal. 
Hereby, the points in time in which the extreme values 
30 are located may be a direct measure of the heart rate. 

This is feasible in particular when the post processor is 
further adapted for validating the quality of the signal 
representing the conformity of the biased auscultation 
35 signal. The quality of the estimated heart rate may 
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hereby be taken into consideration when the estimated 
heart rate is interpreted. 

The quality may be validated by verifying at least one of 
5 the following three items: 

a) time differences between located extreme values must 
be within predetermined limits; 

b) minimum and maximum time differences in proportion to 
the mean of the time differences must be within prede- 

10 termined limits; 

c) the magnitude of the result of the correlation at the 
. extreme values location must be within predetermined 

limits . 

Then a preferred quality level may be defined by means of 
15 the preferred limits that is related directly to the sta- 
tistics of rhythm estimation. 

The bias processor may comprise an envelope detector for 
calculating a first envelope signal representing the en- 

20 velope of the A-weighted version of the auscultation sig- 
nal. This is expedient because the envelope signal of the 
A-weighted version of the auscultation signal comprises 
information that is very feasible for estimating the 
heart rate even when murmurs are present in the ausculta- 

25 tion signal. 

When the bias processor comprises an envelope detector 
for calculating a second envelope signal representing the 
envelope of the auscultation signal, the envelope signal 
30 of the auscultation signal is very feasible for estimat- 
ing the heart rate when murmurs are not present in the 
auscultation signal. 

The bias processor may comprise an adaptive band-pass 
35 filter for filtering the envelope signals; said band-pass 
filter at least having an upper and a lower pass-band. 
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respectively, to select from; said adaptive band-pass 
filter comprising a controller selecting the lower pass- 
band when a relatively large fraction of a signal input 
to the band-pass filter is low-frequent, and selecting 
the upper pass-band when a relatively low fraction of a 
signal input to the band-pass filter is low-frequent. 
Thus, it is possible to estimate the rhythm in the aus- 
cultation signal over a relatively wide range of frequen- 
cies while it is avoided that the signal-to-noise ratio 
is ruined, i.e. it is possible preserve a relatively nar- 
row pass-band of the band-pass filter. 



The bias processor may comprise means for combining the 
envelope signals or the band-pass filtered envelope sig- 
15 nals; said means comprising auto-scaling means and summa- 
tion means providing the biased auscultation signal. It 
is thereby possible to combine differently biased signals 
in order to estimate the heart rate of a broader spectrum 
of types of auscultation signals. 

20 

The auscultation signal may comprise samples that arrive 
at a sample rate and the apparatus may comprise a syn- 
chronous processor operating at a rate corresponding to 
the sample rate, and further comprising an asynchronous 
25 processor operating at time intervals that are initiated 
by a request. Hereby, it is possible to acquire new sam- 
ples of the auscultation signal in real time while a cor- 
relation function operating on prior sample values of the 
auscultation signal is calculated. 

30 

A third object of the invention is to provide a method 
for processing auscultation signals in order to extract 
information representative of the rhythm of a beating 
heart such that random noise and spontaneous artefacts 
35 are suppressed. 
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The quality may be validated by verifying at least one of 
the following three items: 

a) time differences between located extreme values must 
be within predetermined limits; 
5 b) minimum and maximum time differences in proportion to 
the mean of the time differences must be within prede- 
termined limits; 
c) the magnitude of the result of the correlation at the 
extreme values location must be within predetermined 
10 limits. 

Then a preferred quality level may be defined by means of 
the preferred limits that are related directly to the 
statistics of rhythm estimation. 

15 The method may further comprise the step of calculating 
an A-weighted version of the auscultation signal or an 
approximated A-weighted version of the auscultation sig- 
nal. Simulations have shown that an A-weighted version of 
the auscultation signal comprises sufficient information 

20 to estimate the heart rate in a consistent way in a vari- 
ety of situations. 

When the A-weighted version of the auscultation signal is 
calculated by means of an approximation corresponding to 
25 a double differentiation of the auscultation signal, the 
filtering process requires a minimum of computational ef- 
fort. 

The method may comprise the step of calculating a first 
30 envelope signal representing -the envelope of the A- 
weighted version of the auscultation signal. This is ex- 
pedient because the envelope signal of the A-weighted 
version of the auscultation signal comprises information 
that is very feasible for estimating the heart rate even 
35 when murmurs are present in the auscultation signal. 
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When the method further comprises the step of calculating 
a second envelope signal representing the envelope of the 
auscultation signal, the envelope signal of the ausculta- 
5 tion signal is very feasible for estimating the heart 
rate when murmurs are not present in the auscultation 
signal . 

The method may further comprise the step of filtering the 
10 envelope signals by means of a band-pass filter; said 
band-pass filter at least having an upper and a lower 
pass-band respectively selectable; said adaptive band- 
pass filter being controlled such that the lower pass- 
band is selected when a relatively large fraction of a 
15 signal input to the band-pass filter is low-frequent and 
such that the upper pass-band is selected when a rela- 
tively low fraction of a signal input to the band-pass 
filter is low-frequent. Thus, it is possible to estimate 
the rhythm in the auscultation signal over a relatively 
20 wide range of frequencies while it is avoided that the 
signal-to-noise ratio is ruined, i.e. it is possible pre- 
serve a relatively narrow pass-band of the band-pass fil- 
ter . 

25 The method may further comprise the steps of auto-scaling 
the envelope signals and adding the auto-scaled envelope 
signals in order provide the biased auscultation signal. 
It is thereby possible to combine differently biased sig- 
nals in order to estimate the heart rate of a broader 

30 spectrum of types of auscultation signals. 

The auscultation signal may comprise samples that arrive 
at a sample rate and in that the method comprises syn- 
chronous steps being executed at a rate corresponding to 
35 the sample rate, and further comprising asynchronous 
steps operating at time intervals that are initiated by a 
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request. Hereby, it is possible to acquire new samples of 
the auscultation signal in real time while a correlation 
function operating on prior sample values of the auscul- 
tation signal is calculated. 

5 

A fourth object of the invention is to provide a method 
for processing auscultation signals in order to extract 
information relating to the repeated interdependent con- 
traction and expansion of a beating heart, such that the 
10 information in the characteristic sounds and in murmurs 
related to the heart rate are enhanced. 

This is achieved, when the method is characterised in 
that the signal representative of the estimated rhythm is 
15 calculated by means of a filter having a frequency re- 
sponse corresponding to an A-weighing or an approximated 
A-weighing, at least for a frequency range of interest. 

Consequently, it is possible to estimate the heart rate 
20 of a beating heart, based on the characteristic sounds of 
the beating heart and murmurs related to the heart rate. 
Thus, the heart rate may be estimated based on murmurs 
relating to the heart rate if the characteristic sounds 
do not provide sufficient information for estimating the 
25 heart rate. 

When the frequency response is obtained by means of a 
double differentiation, the computational effort required 
to calculate the A-weighted version of the auscultation 
30 is reduced. 

The frequency range of interest may span the frequency 
range below 2000Hz. The characteristic heart sounds and 
murmurs related to the heart rate comprise frequency com- 
35 ponents located in this frequency range and provide suf- 
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ficient information in this range for estimating the 
rhythm of a beating heart. 

In an expedient method for processing auscultation sig- 
5 nals the method comprises the step of selecting a part of 
the biased auscultation -signal and measuring the confor- 
mity between the part of the biased auscultation signal 
and the entire auscultation signal, 

10 In an expedient apparatus for processing auscultation 
signals, the method calculates a cross-correlation func- 
tion. This cross-correlation function supplies a well- 
defined signal having large signal values at good confor- 
mity and small signal values at poor conformity. 

15 

When the part of the biased auscultation signal repre- 
sents one of a succession of cycles of the biased auscul- 
tation signal/ the part may represent a reference cycle 
that is found ideal in some aspects. 

20 

The method may further comprise the step of locating ex- 
treme values in the signal representing the conformity of 
the biased auscultation signal. Hereby, the points in 
time in which the extreme values are located may be a di- 
25 rect measure of the heart rate. 

This is feasible in particular when the method further 
comprises a step of validating the quality of the signal 
representing the conformity of the biased auscultation 
30 signal. The quality of the estimated heart rate may 
hereby be taken into consideration when the estimated 
heart rate is interpreted. 

The quality may be validated by verifying at least one of 
35 the following three items: 
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a) time differences between located extreme values must 
be within predetermined limits; 

b) minimum and maximum time differences in proportion to 
the mean of the time differences must be within prede- 

5 termined limits; 

. c) the magnitude of the result of the correlation at the 
extreme values location must be within predetermined 
limits . 

Then a preferred quality level may be defined by means of 
10 the preferred limits that is related directly to the sta- 
tistics of rhythm estimation. 

The method may comprise the step of calculating a first 
envelope signal representing the envelope of the A- 
15 weighted version of the auscultation signal. This is ex- 
pedient because the envelope signal of the A-weighted 
version of the auscultation signal comprises information 
that is very feasible for estimating the heart rate even 
when murmurs are present in the auscultation signal- 

20 

When the method comprises the step of calculating a sec- 
ond envelope . signal representing the envelope of the aus- 
cultation signal, the envelope signal of the auscultation 
signal is very feasible for estimating the heart rate 
25 when murmurs are not present in the auscultation signal. 

The method may further comprise the step of filtering the 
envelope signals by means of a band-pass filter; said 
band-pass filter at least having an upper and a lower 

30 pass-band respectively selectable; said adaptive band- 
pass filter being controlled such that the lower pass- 
band is selected when a relatively large fraction of a 
signal input to the band-pass filter is low-frequent and 
such that the upper pass-band is selected when a rela- 

35 tively low fraction of a signal input to the band-pass 
filter is low-frequent. Thus, it is possible to estimate 
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the rhythm in the auscultation signal over a relatively 
wide range of frequencies while it is avoided that the 
signal-to-noise ratio is ruined, i.e. it is possible pre- 
serve a relatively narrow pass-band of the band-pass fil- 
5 ter. 

The method may comprise the step of auto-scaling the en- 
velope signals and adding the auto-scaled envelope sig- 
nals in order to provide the biased auscultation signal • 
10 It is thereby possible to combine differently biased sig- 
nals in order to estimate the heart rate of a broader 
spectrum of types of auscultation signals. 

The auscultation signal may comprise samples that arrive 
15 at a sample rate and in that the method comprises syn- 
chronous steps being executed at a rate corresponding to 
the sample rate, and further comprising asynchronous 
steps being executed at time intervals that are initiated 
by a request. Hereby, it is possible to acquire new sam- 
20 pies of the auscultation signal in real time while a cor- 
relation function operating on prior sample values of the 
auscultation signal is calculated. 

Brief Description of the Drawings 
25 The invention will be explained more fully below in con- 
nection with a preferred embodiment and with reference to 
the drawing, in which: 

fig. 1 shows a block diagram of a processor for process- 
30 ing auscultation signals by means of a bias filter calcu- 
lating an approximated A-weighted auscultation signal; 
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fig. 2 shows the expected frequency spectrum of the char- 
acteristic heart sounds, a murmur, and rumble noise; 
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fig. 3 shows a block diagram of a processor for calculat- 
ing the beat frequency of an auscultation signal by means 
of a correlator; . 

5 fig. 4 shows a block diagram of a processor for calculat- 
ing the beat frequency of an auscultation signal by means 
of a correlator and a bias filter comprising A-weighing; 

fig. 5 shows an adaptive band-pass filter; 

fig. 6 shows signals related to an embodiment using a 
correlator and a bias filter comprising A-weighing; 

fig. 7 shows an example of the information in the enve- 
15 lope of the auscultation signal being insufficient; 

fig. 8 shows an example of the information in the enve- 
lope of the A-weighted auscultation signal being insuffi- 
cient. 

20 

Detailed Description of Preferred Embodiments 
Fig. 1 shows a block diagram of a processor for process- 
ing auscultation signals by means of a bias filter calcu- 
lating an approximated A-weighted auscultation signal. An 

25 auscultation signal is a signal representing the acousti- 
cal signal that may be picked up by some type of trans- 
ducer located to provide an auscultation signal repre- 
senting the sound of a beating heart. The auscultation 
signal is hereinafter denoted VCG and may be a digital or 

30 analogue signal. 

The block diagram comprises a bias processor 106. The VCG 
signal is filtered in the bias filter 101. The frequency 
response of the filter corresponds to an approximated A- 
35 weighing of the VCG signal. A-weighing is a linear model 
of the human ear modelling the human perception of signal 
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amplitudes in the audio frequency range, i.e. about 20- 
20000Hz, A-weighing is typically modelled as a second or- 
der band-pass filter with a centre frequency of 2000Hz. 
The A-weighing model is approximated as a double differ- 
5 entiation of the VCG signal with respect to time thereby 
providing an A-weighted version, A-VCG, of the VCG sig- 
nal. The approximation is valid for frequencies between 
about 20-200H2. The bias filter 101 may thus be viewed as 
a second order high-pass filter slope. The filter thereby 

10 diminishes low frequency rumble noise that may be present 
in the auscultation signal. Further, the filter enhances 
the upper part of the frequency spectrum of the charac- 
teristic heart sounds and eventually murmurs. However, 
the A-weighted signal A-VCG may also enhance noise origi- 

15 nating from sources that are not related to either the 
characteristic heart sounds or eventually murmurs thereby 
destroying the signal-to-noise-ratio. The A-weighted sig- 
nal A-VCG is therefore low-pass filtered by means of the 
filter 102. The low-pass filtered signal is then supplied 

20 as an input signal to the rectifier 102 that calculates 
the absolute value of the input signal and thereby pro- 
vides an output signal. The output signal from the recti- 
fier 103 is band-pass filtered in the filter 104, This 
band-pass filtered signal thereby constitutes the enve- 

25 lope of the A-weighted signal A-VCG. The band-pass filter 
has a pass-band corresponding to the range of expected 
heart rates e.g, from 0,5 to 3,3HZ, Further, the stop- 
band of the band-pass filter diminishes any DC component 
and frequency components above the fundamental heart beat 

30 frequency, i,e. relative high frequency components. 

The information in the band-pass filtered signal present 
as repetitive well-defined cycles is thereby representa- 
tive of the heart rate. The heart rate may be estimated 
35 by means of an estimator 105. The estimator may be in the 
form of a peak locater determining the timing of the 
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well-defined extreme values of the band-pass filtered 
signal. Alternatively, the estimator may be in the form 
of a threshold or zero-crossing detector thereby generat- 
ing a square signal with a frequency representing the es- 
5 timated heart rate. When the estimator comprises a 
threshold detector, the biased auscultation signal and a 
threshold value are both input to a comparator. The 
threshold value may be in the form of a constant value or 
an adaptive threshold value. The adaptive threshold _ value 
10 may be provided by a low-pass filter filtering the biased 
auscultation signal . 

Fig. 2 shows the expected frequency spectrum of the char- 
acteristic heart sounds, a murmur, and rumble noise. The 

15 spectrums are rough approximations and for illustration 
purposes only. The magnitude Amp. of the spectrums are in 
arbitrary units and shown as a function of the frequency 
F [Hz] . The frequency spectrum 301 of the characteristic 
heart sounds related to a contraction and an expansion of 

20 a beating heart is illustrated roughly. In the literature 
the sound related to a contraction is denoted SI and the 
sound related to an expansion is denoted S2. The fre- 
quency spectrum of the frequency contents of SI and S2 is 
primarily located between 20 and 50Hz but may be distrib- 

25 uted up to 100-200HZ and down to below lOHz. 

However, the frequency contents of murmurs are typically 
distributed at higher frequencies. This is shown by the 
frequency response 302. The frequency contents of murmurs 

30 comprising ^musical' sounds is distributed over a rela- 
tively narrow-band frequency range 306. This relative 
narrow band comprises information that may be used for 
estimation of the heart rate. Typically, murmurs are also 
associated with a relative broad band 205 not comprising 

35 information regarding the heart rate. An absolute numeri- 
cal definition of these frequency bands and the magni- 
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tudes of the auscultation signal in these bands will de- 
pend on the murmur in question . 

The frequency spectrum 304 of low-frequent rumble noise 
5 typically overlaps the frequency spectrum 301 of the 
. characteristic heart sounds. 

The transfer function 303 shows a preferred transfer 
function for a filter for calculating a biased ausculta- 

10 tion signal. The filter is thereby able to enhance the 
amplitude of the frequency components of the upper part 
of the frequency spectrum of the characteristic heart 
sounds and the amplitude of the frequency components of 
the frequency spectrum of murmurs. Further, low-frequency 

15 rumble noise that may be present in the auscultation sig- 
nal is diminished. The transfer function 303 comprises a 
high-pass 40 dB per decade slope. The slope may be a part 
of an approximated A-weighing function, where the ap- 
proximation is implemented as a second order band-pass 

20 filter with a centre frequency of 2000Hz. 

Fig. 3 shows a block diagram of a processor for calculat- 
ing the beat frequency of an auscultation signal by means 
of a correlator. In this embodiment the received auscul- 
25 tation signal VCG is a digital signal comprising digital 
samples arriving at a sample frequency fs, e.g. fs equal 
to 4000Hz. 

The embodiment comprises a bias processor 201 processing 
30 the VCG samples at a rate synchronous with the sample 
frequency and an estimator 302 capable of operating asyn- 
chronous relative to the bias processor and providing a 
signal HR representative of the heart rate, when it is 
requested- The estimator 302 receives data in the form of 
35 samples from the bias processor 301 via the buffer 303. 
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The buffer 303 may be in the form of a circular buffer or 
a first-in-first-out (FIFO) memory. 

In the bias processor 301 the digital auscultation signal 
5 VCG is received and the respective absolute values of the 
VCG signal are calculated by means of the rectifier 304. 
Then the absolute values are band-pass filtered by means 
of the filter 305 in order to calculate the envelope of 
the auscultation signal. The band-pass filter has a _pass- 
10 band corresponding to the range of expected heart rates 
e.g. from 0.5 to 3 . 3HZ . Further, the stop-band of the 
band-pass filter diminishes any DC component and fre- 
quency components above the fundamental heart beat fre- 
quency, i.e. relatively high frequency components. 

15 

Samples from the band-pass filter are supplied to the 
buffer 303 at even points in time. The buffers may e.g. 
be circular buffers that have a capacity for a number of 
sample points, e.g corresponding to a signal sequence of 

20 5 seconds. The oldest sample in the buffer is overwritten 
when a new sample arrives from the band-pass filter. Sam- 
ples are taken from the buffer when the estimator 302 is 
initiated. The estimator may be initiated at certain 
regular or irregular points in time e.g.. in order to pro- 

25 vide information for the update of a digital display. 

Although the envelope signal provided by the band-pass 
filter via the buffer 303 is smooth and comprises repeti- 
tive cycles it may be difficult to determine whether one 

30 or more cycles are present during a specific time inter- 
val. However, the information content in the envelope 
signal may be enhanced towards determining the exact num- 
ber of cycles in a specific time interval. This is done 
by calculating the conformity of the envelope signal with 

35 respect to the envelope signal itself. This is based on 
the assumption that the vast majority of the cycles in 



25 

the envelope signal represents an actual interrelated 
contraction and expansion of a heart in the process of 
pumping blood. The conformity of the envelope signal is 
calculated by means of the correlator 306 calculating the 
5 auto-correlation of the envelope signal. 

In a preferred embodiment the correlator comprises means 
for selecting that cycle in the envelope signal that has 
the largest amplitude. This is done under the fair as- 

10 sumption that the largest amplitude in the auscultation 
signal is caused by a normal contraction/expansion of the 
heart that is least likely to comprise any artefacts. The 
extent of the selected cycle is found by searching for 
the nearest local minima that is less than zero. Thereby, 

15 falling into, local minima with a magnitude greater than 
zero is avoided. This search for a minima is applied both 
forward and backward relative to the position of the 
largest amplitude. The selected cycle is stored and 
cross-correlated with the rest of the envelope signal. 

20 Alternatively, convolution or any other function capable 
calculating the conformity of a signal with respect to 
the signal itself may be used. 

Fig. 4 shows a block diagram of a processor for calculat- 
25 ing the beat frequency of an auscultation signal by means 
of a correlator and a bias filter comprising A-weighing. 
In this preferred embodiment the bias filter comprising 
A-weighing and the correlator are combined in order to 
obtain an enhanced overall performance, i.e. it is possi- 
30 ble to estimate the rhythm of a beating heart even though 
the characteristic sounds of the heart are distorted se- 
verely. 

The auscultation signal VCG is in the form of a digital 
35 signal comprising sample values arriving at a sample fre- 
quency fs, e.g. fs = 4000Hz. The auscultation signal is 
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branched. A first branch supplies the auscultation signal 
to a rectifier or an absolute value detector, ABS, 406. 
The absolute value of the auscultation signal is further 
supplied to the low-pass filter 407 having a cut-off fre- 
5 quency of about 3Hz such that frequency components above 
lOOHz is diminished. A second branch supplies the auscul- 
tation signal to the filter 401 in which the double de- 
rivative of the auscultation signal is calculated. The 
frequency response of the filter 401 corresponds . to a 

10 second order high-pass slope or an A-weighing approxima- 
tion in the frequency range below 2000Hz. The double de- 
rivative of the auscultation signal is supplied to a low- 
pass filter 402 in which the relatively high frequency 
noise components enhanced by the differentiation in the 

15 filter 401 is damped. The filter 402 may have a cut-off 
frequency of about lOOHz. The signal that is low-pass 
filtered by the filter 402 is supplied to a rectifier or 
an absolute value detector, ABS, 403. The absolute value 
of the auscultation signal is further supplied to the 

20 low-pass filter 407 that diminishes frequency components 
above lOOHz. 

The decimators, DEC, 405 and 408 generate a reduced sam- 
ple rate fsa, that e.g. may be 200Hz since frequency com- 
25 ponents above lOOHz is diminished. Thus, the computa- 
tional effort required in the succeeding signal process- 
ing of the auscultation signal may be reduced. 

The band-pass filter 411 comprising the band-pass filters 
30 BPF#1 and BPF#2 receives the decimated auscultation sig- 
nals and calculates band-pass filtered signals. The band- 
pass filter may have a pass-band corresponding to the 
range of expected heart rates e.g. from 0.5 to 3.3HZ. 
Note that DC components and frequency components above 
35 25Hz are diminished. The centre frequency of the band- 
pass filter may be adjusted adaptively in order to en- 



27 



hance the performance of the rhythm estimator. This is 
discussed later. 

The band-pass filtered signals thus represent the enve- 
5 lopes of the auscultation signal and the double differen- 
tiated auscultation signal at a reduced sample rate. 

The decimators, DEC, 412 and 413 generate a further re- 
duced sample rate that may be e.g. 50Hz since frequency 
10 components above 25Hz are diminished. Thus, the computa- 
tional effort required in the succeeding signal process- 
ing of the auscultation signal may be reduced even fur- 
ther . 

15 Samples from the band-pass filter are supplied to the 
buffers 414 and 415 at even points in time. The buffers 
may e.g. be circular buffers that have a capacity for a 
number of sample points, e.g corresponding to a signal 
sequence of 5 seconds. Samples are taken from the buffers 

20 when the estimator 422 is initiated. The estimator may be 
initiated at certain regular or irregular points in time 
e.g. in order to provide information for the update of a 
digital display. 

25 The samples that are taken from the buffers 414 and 415 
are scaled by means of the auto-scale means 416 and 417 
respectively. The auto-scale means may scale the received 
sample values by ensuring that the peak value over a num- 
ber of scaled samples take the value 1.0, e.g. the number 

30 of samples may correspond to the signal sequence stored 
in the buffers 414 and 415, that is the signal sequence 
is normalised to a peak value equal to 1.0. The scaled 
samples are added in the adder 418 with equal weights 
since the information in the two signals are of equal im- 

35 portance, and the result of the addition is provided as a 
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sum-signal. Alternatively, the two signals may be added 
with different weights. 

The conformity of the sum-signal is calculated by means 
5 of the correlator, X-Corr, 419. In a preferred embodiment 
the correlator comprises means for selecting that cycle 
in the envelope signal that has the largest amplitude. 
The time extension of the selected cycle is found by 
searching for the nearest local minima that is less, than 

10 zero. Thereby, falling into local minima with a magnitude 
greater than zero is avoided. This search for a minima is 
applied both forward and backward relative to the posi- 
tion of the largest amplitude. The selected cycle is 
stored and cross-correlated with the rest of the envelope 

15 signal. 

The peak-locater 420 identifies extreme values in a sig- 
nal representative of the result of the cross- 
correlation. This may e.g. be carried out by calculating 
20 the first derivative of the result of the cross- 
correlation and identifying inflexions. Preferably the 
peak-locater identifies maximums in the signal represen- 
tative of the result of the cross-correlation. 

25 The locations of the extreme values are preferably proc- 
essed in a quality estimator 421. The estimator decides 
whether the result of the cross-correlation qualifies for 
estimating the rhythm, if any, of the auscultation sig- 
nal. The quality estimator 421 may comprise three steps 

30 in the process of estimating the quality: 

1. The mean of the time differences between two extreme 
values must be within predetermined limits; 

2. The minimum and maximum time differences in proportion 
to the mean of the time differences must be within 

35 predetermined limits; 
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3. The magnitude of the result of the cross-correlation 
at the extreme values location must be within prede- 
termined limits. 



5 If one or more of the above mentioned items 1-3 fail, 
. then it is determined that, it is not possible to estimate 
any rhythm in the auscultation signal. Alternatively, if 
the above mentioned items pass, then it is possible to 
. estimate a rhythm in the auscultation signal. The rhythm 
10 may be estimated by counting the number of extreme values 
in the result of the cross-correlation per time unit and 
displaying this as a number of beats per minute . 

Fig. 5 shows an adaptive band-pass filter. This filter 
15 may be used instead of the band-pass filter 411 shown in 
fig. 4. The adaptive filter may be used in estimation of 
the rhythm in the auscultation signal over a relatively 
wide range of frequencies while it is avoided that the 
signal-to-noise ratio is destroyed, i.e. it is possible 
20 preserve a relatively narrow pass-band of the band-pass 
filter. 

The band-pass filter 512 receives two input signals II 
and 12 and provides two band-pass filtered output signals 

25 01 and 02 . The band-pass filter 512 thus comprises two 
band-pass filters BPF#1 and BPF#2 . Each of the filters 
BPF#1 and BPF#2 comprises a pass-band LB with a rela- 
tively low centre frequency and a pass-band HB with a 
relatively high centre frequency. One of the pass-bands 

30 may be selected by supplying an enable signal to either 
the input LB or HB . This pass-band selection is carried 
out by the adaptive controller 513, 
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The high-pass filter 501 receives an input signal corre- 
sponding to one of the signals input to the band-pass 
filter 512. The high-pass filter 501 has a cut-off fre- 
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quency of about 0.3Hz thereby diminishing very low fre- 
quency components that is not used in the adaptive pass- 
band selection. The filter 501 supplies a high-pass fil- 
tered signal to the low-pass filter 502. The low-pass 
5 filter has a cut-off frequency of about IHz . 

The buffers 503 and 504 are capable of storing sample 
values corresponding to a signal sequence with a time du- 
ration of about 3 seconds, respectively. The expected or 

10 mean value of the sequences stored in the buffers are 
calculated by means of the mean value detectors 505 and 
506. The mean value calculated by mean value detector 505 
represents the mean value of the entire signal input to 
the band-pass filter 512. The mean value calculated by 

15 mean value detector 506 represents the mean value of the 
low-frequency part, i.e. below about IHz, of the signal 
input to the band-pass filter 512. 

Thus, the. divider 507 calculates the fraction A/B of the 
20 signal that is relatively low-frequent. In an initial 
situation the threshold Th. is equal to the value 0.45.. 
This value is subtracted from the fraction A/B by means 
of the adder 508. In the decision state 509 it is veri- 
fied whether the frequency ratio A/B is greater than the 
25 threshold value. If this is the case (T) then the thresh- 
old value Th, is set equal to the value 0.35 and the 
pass-band LB having a relatively low centre frequency is 
selected. This is carried out under the assumption that 
the auscultation signal originates from a slow beating 
30 heart. Thus it is possible to diminish higher and super- 
fluous frequency components by selecting the pass-band LB 
having a relatively low centre frequency. Preferably, the 
relatively low centre frequency is about 0,7Hz. 

35 Alternatively, if the frequency ratio A/B is less than 
the threshold value Th.(F), then the threshold value Th. 
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is set equal to the value 0.55 and the pass-band HB hav- 
ing a relatively high centre frequency is selected. This 
is carried out under the assumption that due to the lack 
of low-frequency content in the auscultation signal, the 
5 auscultation is likely to originate from a very fast 
beating heart, e.g. from, a child. By selecting the pass- 
band HB having a relatively high centre frequency there 
is maximum likelihood that the rhythm of the auscultation 
signal is estimated correctly. Preferably, the relatively 
10 high centre frequency is about 1.2Hz. 

Fig. 6 shows signals related to an embodiment using a 
correlator and a bias filter comprising A-weighing. The 
auscultation signal VCG and the corresponding band-pass 

15 filtered envelope of the VCG signal is shown as a func- 
tion of time. Further, the A-weighted or double differen- 
tiated auscultation signal d^VCG/dt^ and the correspond- 
ing band-pass filtered envelope is shown. In view of a 
graphical interpretation it is clear that the heart rate 

20 could be estimated by using only one of the band-pass 
filtered signals. However, a stronger estimate is pro- 
vided if the signals are added and used as input to a 
correlator as described above. This is shown by means of 
the signals 'XC-I' and ^XC-0' . 

25 

Fig. 7 shows an example of when the information in the 
envelope of the auscultation signal is insufficient. The 
auscultation signal VCG and the corresponding band-pass 
filtered envelope of the VCG signal is shown as a func- 

30 tion of time. Further, the A-weighted or double differen- 
tiated auscultation signal drVCG/dt^ and the correspond- 
ing band-pass filtered envelope is shown. In view of a 
graphical interpretation it is clear that the strongest 
estimate is provided by the double differentiated auscul- 

35 tation signal d^VCG/dt^. -However, for the vast majority 



of auscultation signals the signal 'XC-O' provides the 
strongest overall estimate. 

Fig. 8 shows an example of when the information in the 
envelope of the A-weighted auscultation signal is insuf- 
ficient. The auscultation signal VCG and the correspond-- 
ing band-pass filtered envelope of the VCG signal, is 
shown as a function of time. Further, the A-weighted or 
double differentiated auscultation signal d^VCG/dtf and 
the corresponding band-pass filtered envelope is shown. 
In view of a graphical interpretation it is clear that 
the band-pass filtered VCG envelope signal provides the 
strongest estimate. But, again, the signal 'XC-0' pro- 
vides the strongest overall estimate. 

The invention may be embodied as a part of a stethoscope 
or any other instrument or apparatus. Alternatively, the 
invention may be embodied as a computer program or a part 
of a computer program, which may be loaded into the mem- 
ory of a computer and executed therefrom. The computer 
program may be distributed by means of any data storage 
or data transmission medium. The storage media can be 
magnetic tape, optical disc, compact disc (CD or CD-ROM), 
mini-disc, hard disk, floppy disk, ferroelectric memory, 
electrically erasable programmable read only memory 
(EEPROM) , flash memory, EPROM, read only memory (ROM) , 
static random access memory (SRAM) , dynamic random access 
memory (DRAM) , ferromagnetic memory, optical storage, 
charge coupled devices, smart cards, etc. The transmis- 
sion medium can be a network, e.g. a local area network 
(LAN) , a wide area network (WAN) , or any combination 
thereof, e.g. the Internet. The network may comprise wire 
and wire-less communication links. Via the network a 
software embodiment (i.e. a program) of the invention, or 
a part thereof, may be distributed by transferring a pro- 
gram via the network. 
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CLAIMS 

1. An apparatus for processing auscultation signals, com- 
prising 

5 

a bias processor (106; 301; 409) for receiving an auscul- 
tation signal and providing a biased auscultation signal; 
said bias processor comprising an envelope detector; 

10 an estimator (105; 302; 410) for calculating a signal 
representative of an estimated rhythm of the auscultation 
signal; 

characterised in that, 

15 

the estimator is adapted for selecting a part of the bi- 
ased auscultation signal and measuring the conformity be- 
tween the part of the biased auscultation signal and the 
entire auscultation signal. 

20 

2. An apparatus for processing auscultation signals ac- 
cording to claim 1, characterised in that the estimator 
calculates a cross-correlation function, 

25 3. An apparatus for processing auscultation signals ac- 
cording to claim 1-2, characterised in that the part of 
the biased auscultation signal represents one of a suc- 
cession of cycles of the biased auscultation signal, 

30 4 . An apparatus for processing auscultation signals ac- 
cording to claim 1-3, characterised in that the apparatus 
further comprises a post processor (420, 421) adapted for 
locating extreme values in the signal representing the 
conformity of the biased auscultation signal. 
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5. An apparatus for processing auscultation signals ac- 
cording to claim 4, characterised in that the post proc- 
essor (420, 421) further is adapted for validating the 
quality of the signal representing the conformity of the 
biased auscultation signal. 

6. An apparatus for processing auscultation signals ac- 
cording to claim 5, characterised in that the quality is 
validated by verifying at least of the following ihree 
items: 

a) time differences between located extreme values must 
be within predetermined limits; 

b) minimum and maximum time differences in proportion to 
the mean of the time differences must be within prede- 
termined limits; 

c) the magnitude of the result of the correlation at the 
extreme values location must be within predetermined 
limits. 

7. An apparatus for processing auscultation signals ac- 
cording to claim 1-6, characterised in that the bias 
processor comprises a filter (101; 401) for calculating 
an A-weighted version of the auscultation signal or an 
approximated A-weighted version of the auscultation sig- 
nal . 

8. An apparatus for processing auscultation signals ac- 
cording to claim 7, characterised in that the A-weighted 
version of the auscultation signal is calculated by means 
of an approximation corresponding to a double differen- 
tiation of the auscultation signal. 

9. An apparatus for processing auscultation signals ac- 
cording to claim 1-8, characterised in that the bias 
processor comprises an envelope detector (102,103; 
403,404) for calculating a first envelope signal repre- 



senting the envelope of the A-weighted version of the 
auscultation signal • 

10. An apparatus for processing auscultation signals ac- 
cording to claim 1-9, characterised in that the bias 
processor (302; 409) comprises an envelope detector 
(304,305; 406,407) for calculating a second envelope sig- 
nal representing the envelope of the auscultation signal. 

11. An apparatus for processing auscultation signals ac- 
cording to claim 9-10, characterised in that the bias 
processor comprises an adaptive band-pass filter (512) 
for filtering the envelope signals; said band-pass filter 
at least having an upper and a lower pass-band respec- 
tively selectable; said adaptive band-pass filter com- 
prising a controller (513) selecting the lower pass-band 
when a relatively large fraction of a signal input to the 
band-pass filter is low-frequent and selecting the upper 
pass-band when a relatively low fraction of a signal in- 
put to the band-pass filter is low-frequent. 

12. An apparatus for processing auscultation signals ac- 
cording to claim 9-11, characterised in that the bias 
processor comprises means for combining the envelope sig- 
nals or the band-pass filtered envelope signals; said 
means comprising auto-scaling means (416, 417) and summa- 
tion means (418) providing the biased auscultation sig- 
nal . 

13. An apparatus for processing auscultation signals ac- 
cording to claim 1-12, characterised in that the auscul- 
tation signal comprises samples that arrive at a sample 
rate and in that the apparatus comprises a synchronous 
processor (301; 409) operating at a rate corresponding to 
the sample rate, and further comprising an asynchronous 



processor (302; 410) operating at time intervals that are 
initiated by a request. 

14. An apparatus for estimating the rhythm in ausculta- 
5 tion signals, comprising 

a bias processor (106; 301; 409) for receiving an auscul- 
tation signal and providing a biased auscultation signal; 
said bias processor comprising an envelope detector;. 

10 

an estimator (105; 302; 410) for calculating a signal 
representative of an estimated rhythm of the auscultation 
signal; 

15 characterised in that, 

the bias processor comprises a filter (101; 401) having a 
frequency response corresponding to an A-weighing or an 
approximated A-weighing, at least for a frequency range 
20 of interest- 

15. An apparatus for processing auscultation signals ac- 
cording to claim 14, characterised in that the ausculta- 
tion signal is filtered with a filter (101; 401) having a 

25 frequency response corresponding to a double differentia- 
tion. 

16. An apparatus for processing auscultation signals ac- 
cording to claim 14-15, characterised in that the fre- 

30 quency range of interest spans the frequency range below 
2000Hz. 

17. An apparatus for processing auscultation signals ac- 
cording to claim 14-16, characterised in that the estima- 

35 tor is adapted for selecting a part of the biased auscul- 
tation signal and measuring the conformity between the 



part of the biased auscultation signal and the entire 
auscultation signal. 



18. An apparatus for processing auscultation signals ac- 
cording to claim 17, characterised in that the estimator 
calculates a cross-correlation function. 

19. An apparatus for processing auscultation signals ac- 
cording to claim 17-18, characterised in that the part of 
the biased auscultation signal represents one of a suc- 
cession of cycles of the biased auscultation signal. 

20. An apparatus for processing auscultation signals ac- 
cording to claim 14-19, characterised in that the appara- 
tus further comprises a post processor (420, 421) adapted 
for locating extreme values in the signal representing 
the conformity of the biased auscultation signal. 

21. An apparatus for processing auscultation signals ac- 
cording to claim 20, characterised in -that the post proc- 
essor (420, 421) further is adapted for validating the 
quality of the signal representing the conformity of the 
biased auscultation signal. 

22. An apparatus for processing auscultation signals ac- 
cording to claim 21, characterised in that the quality is 
validated by verifying at least one of the following 
three items: 

a) time differences between located extreme values must 
be within predetermined limits; 

b) minimum and maximum time differences in proportion to 
the mean of the time differences must be within prede- 
termined limits; 

c) the magnitude of the result of the correlation at the 
extreme values location must be within predetermined 
limits . 
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23. An apparatus for processing auscultation signals ac- 
cording to claim 14-22, characterised in that the bias 
processor comprises an envelope detector (102,103; 
403,404) for calculating a first envelope signal repre- 
senting the envelope of the A-weighted version of the 
auscultation signal . 

24. An apparatus for processing auscultation signals ac- 
cording to claim 14-23, characterised in that the bias 
processor comprises an envelope detector (304,305;. 
406,407) for calculating a second envelope signal repre- 
senting the envelope of the auscultation signal. 

25. An apparatus for processing auscultation signals ac- 
cording to claim 23-24, characterised in that the bias 
processor, comprises an adaptive band-pass filter (512) 
for filtering the envelope signals; said band-pass filter 
at least having an upper and a lower pass-band respec- 
tively selectable; said adaptive band-pass filter com- 
prising a controller (513) selecting the lower pass-band 
when a relatively large fraction of a signal input to the 
band-pass filter is low-frequent and selecting the upper 
pass-band when a relatively low fraction of a signal in- 
put to the band-pass filter is low-frequent. 

26. An apparatus for processing auscultation signals ac- 
cording to claim 24-25, characterised in that the bias 
processor comprises means for combining the envelope sig- 
nals or the band-pass filtered envelope signals; said 
means comprising auto-scaling means (416, 417) and summa- 
tion means (418) providing the biased auscultation sig- 
nal . 

27. An apparatus for processing auscultation signals ac- 
cording to claim 14-26, characterised in that the auscul- 
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tation signal comprises samples that arrive at a sample 
rate and in that the apparatus comprises a synchronous 
processor (301; 409) operating at a rate corresponding to 
the sample rate, and further comprising an asynchronous 
5 processor (302; 410) operating at time intervals that are 
initiated by a request. 

28. A method for processing auscultation signals, com- 
prising the steps of 

10 

receiving an auscultation signal and providing a biased 
auscultation signal; 

calculating a signal representative of an estimated 
15 rhythm of the auscultation signal; 

characterised in that, 

the estimated rhythm is calculated by selecting a part of 
20 the biased auscultation signal and measuring the confor- 
mity between the part of the biased auscultation signal 
and the entire auscultation signal. 

29. A method for processing auscultation signals accord- 
25 ing to claim 28, characterised in that the estimator cal- 
culates a cross-correlation function. 

30. A method for processing auscultation signals accord- 
ing to claim 28-29, characterised in that the part of the 

30 biased auscultation signal represents one of a succession 
of cycles of the biased auscultation signal. 

31. A method for processing auscultation signals accord- 
ing to claim 28-30, characterised in that the method fur- 

35 ther comprises the step of locating extreme values in the 
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signal representing the conformity of the biased auscul- 
tation signal. 

32. A method for processing auscultation signals accord- 
5 ing to claim 31, characterised in that the method further 
comprises the step of validating the quality of the sig- 
nal representing the conformity of the biased ausculta- 
tion signal. 

10 33. A method for processing auscultation signals accord- 
ing to claim 32, characterised in that the quality is 
validated by verifying at least one of the following 
three items: 

a) time differences between located extreme values must 
15 be within predetermined limits; 

b) minimum and maximum time differences in proportion to 
the mean of the time differences must be within prede- 
termined limits; 

c) the magnitude of the result of the correlation at the 
20 extreme values location must be within predetermined 

limits . 

34. A method for processing auscultation signals accord- 
ing to claim 28-33, characterised in that the method fur- 

25 ther comprises a step of calculating an A-weighted ver- 
sion of the auscultation signal or an approximated A- 
weighted version of the auscultation signal. 

35. A method for processing auscultation signals accord- 
30 ing to claim 34, characterised in that the A-weighted 

version of the auscultation signal is calculated by means 
of an approximation corresponding to a double differen- 
tiation of the auscultation signal. 



35 36. A method for processing auscultation signals accord- 
ing to claim 28-35, characterised in that the method fur- 



ther comprises the step of calculating a first envelope 
signal representing the envelope of the A-weighted ver- 
sion of the auscultation signal. 

37. A method for processing auscultation signals accord- 
ing to claim 28-36, characterised in that the method fur- 
ther comprises the step of calculating a second envelope 
signal representing the envelope of the auscultation sig- 
nal. 

38. A method for processing auscultation signals accord- 
ing to claim 36-37, characterised in that the method fur- 
ther comprises the step of filtering the envelope signals 
by means of a band-pass filter; said band-pass filter at 
least having an upper and a lower pass-band respectively 
selectable; said adaptive band-pass filter being con- 
trolled such that the lower pass-band is selected when a 
relatively large fraction of a signal input to the band- 
pass filter is low-frequent and such that the upper pass- 
band is selected when a relatively low fraction of a sig- 
nal input to the band-pass filter is low-frequent. 

39. A method for processing auscultation signals accord- 
ing to claim 36-38, characterised in that the method fur- 
ther comprises the steps of auto-scaling the envelope 
signals and adding the auto-scaled envelope signals in 
order provide the biased auscultation signal. 

40. A method for processing auscultation signals accord- 
ing to claim 28-39, characterised in that the ausculta- 
tion signal comprises samples that arrive at a sample 
rate and in that the method comprises synchronous steps 
being executed at a rate corresponding to the sample 
rate, and further comprising asynchronous steps operating 
at time intervals that are initiated by a request. 
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41. A method for estimating the rhythm in auscultation 
signals, comprising the steps of 

receiving an auscultation signal and providing a biased 
auscultation signals- 
calculating a signal representative of an estimated 
rhythm of the auscultation signals- 
characterised in that, 

the signal representative of the estimated rhythm is cal- 
culated by means of a filter having a frequency response 
corresponding to an A-weighing or an approximated A- 
weighing, at least for a frequency range of interest. 

42. A method for processing auscultation signals accord- 
ing to claim 41, characterised in that the frequency re- 
sponse is obtained by means of a double differentiation, 

43. A method for processing auscultation signals accord- 
ing to claim 41-42, characterised in that the frequency 
range of interest spans the frequency range below 2000Hz. 

44. A method for processing auscultation signals accord- 
ing to claim 41-43, characterised in that the method com- 
prises the step of selecting a part of the biased auscul- 
tation signal and measuring the conformity between the 
part of the biased auscultation signal and the entire 
auscultation signal. 

45. A method for processing auscultation signals accord- 
ing to claim 44, characterised in that the method calcu- 
lates a cross-correlation function. 
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46. A method for processing auscultation signals accord- 
ing to claim 44-45, characterised in that the part of the 
biased auscultation signal represents one of a succession 
of cycles of the biased auscultation signal. 

47. ^ method for processing auscultation signals accord- 
ing to claim 41-4 6, characterised in that the. method fur- 
ther comprises a step of locating extreme values in the 
signal representing the conformity of the biased auscul- 
tation signal. 

48. A method for processing auscultation signals accord- 
ing to claim 47, characterised in that the method further 
comprises a step of validating the quality of the signal 
representing the conformity of the biased auscultation 
signal . 

49. A method, for. processing auscultation signals accord- 
ing to claim 48, characterised in that the quality is 
validated by verifying at least one of the following 
three items : 

a) time differences between located extreme values must 
be within predetermined limits; 

b) minimum and maximum time differences in proportion to 
the mean of the time differences must be within prede- 
termined limits; 

c) the magnitude of the result of the correlation at the 
extreme values location must be within predetermined 
limits , 

50. A method for processing auscultation signals accord- 
ing to claim 41-49, characterised in that the method com- 
prises the step of calculating a first envelope signal 
representing the envelope of the A-weighted version of 
the auscultation signal. 
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51, A method for processing auscultation signals accord- 
ing to claim 41-50, characterised in that the method com- 
prises the step of calculating a second envelope signal 
representing the envelope of the auscultation signal. 

52. A method for processing auscultation signals accord- 
ing to claim 50-51, characterised in that the method fur- 
ther comprises the step of filtering the envelope signals 
by means of a band-pass filter; said band-pass filter at 
least having an upper and a lower pass-band respectively 
selectable; said adaptive band-pass filter being con- 
trolled such that the lower pass-band is selected when a 
relatively large fraction of a signal input to the band- 
pass filter is low-frequent and such that the upper pass- 
band is selected when a relatively low fraction of a sig- 
nal input to the band-pass filter is low-frequent. 

53, A method for processing auscultation signals accord- 
ing to claim 51-52, characterised in that the method com- 
prises the step of auto-scaling the envelope signals and 
adding the auto-scaled envelope signals in order to pro- 
vide the biased auscultation signal. 

54. A method for processing auscultation signals accord- 
ing to claim 41-53, characterised in that the ausculta- 
tion signal comprises samples that arrive at a sample 
rate and in that the method comprises synchronous steps 
being executed at a rate corresponding to the sample 
rate, and further comprising asynchronous steps being 
executed at time intervals that are initiated by a re- 
quest . 
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ABSTRACT 



An apparatus for processing auscultation signals, e.g. 
embedded in . an electronic stethoscope with a digital 
readout of an estimated heart rate. The apparatus com- 
prises a bias processor for receiving an auscultation 
signal and providing a biased auscultation signal; said 
bias processor comprising an envelope detector and an es- 
timator for calculating a signal representative of the 
beat frequency of the auscultation signal. In one aspect 
the information in the biased auscultation signal that is 
in conformity with the repeated information in the aus- 
cultation signal is enhanced. This may e.g- be done by 
calculating a cross-correlation between the biased aus- 
cultation signal and a part of the biased auscultation 
signal. In another aspect of the invention the ausculta- 
tion signal is biased by means of A-weighing. Preferably 
the two aspects of the invention are combined. 



Fig. 4 should be published. 



